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Abstract. The variational approach to weakly nonlocal thermodynamic the- 
ories is critically revisited in the light of modern nonequilibrium thermody- 
namics. The example of Ginzburg-Landau equation is investigated in detail. 



1. Introduction 

In the last decades there has been a continuous interest in developing generalized 
classical continuum theories that are able to describe nonlocal effects. One of the 
most important (and popular) strategies is to develop so called weakly nonlocal 
continuum theories that is to incorporate higher order space derivatives into the 
governing equations of continuum physics. The crucial point in this investigations is 
to clarify the relation of the new equations to the Second Law of thermodynamics. 

From that point of view the weakly nonlocal continuum theories can be divided 
into two groups. The theories in the first group take seriously the thermodynamic 
requirements and the estabilished structure of continuum physics. We can call them 
thermodynamic weakly nonlocal continuum theories. The theories in the second 
group does not follow the structure of classical continuum physics and we can call 
them variational weakly nonlocal continuum theories according to the basic method 
of equation construction. Examples of thermodynamic theories are the so called 
gradient theory (of thcrmomcchanics) developed by Kosihski and Valanis [0, ||, ||, 
the virtual power considerations of Germain and Maugin Q , the multifield theory of 
Mariano g[ , a microstructure theory based on the concept of microforcc balance of 
Gurtin [|6| and the investigations toward the weakly nonlocal extension of extended 
thermodynamics by Lebon, Jou and coworkers 0, ^, ^, 0. Instead of reviewing 
and criticizing the different approaches our general remark is, that they usually 
introduce new and disputable new concepts that seem to be too special to serve as 
foundation of a general nonlocal thermodynamic theory. 

The second group is started from the investigations and method introduced by 
Ginzburg and Landau and constructs a set of prototypical classical weakly nonlocal 
equations like the Ginzburg-Landau equation and the Cahn-Hilliard equation (see 
e.g. [0 and the references therein). A unified treatment of weak nonlocality based 
on this variational approach appears in the so called GENERIC scheme developed 
by Grmela and Ottingcr [|2[ ||. 

The following table compares the different nonlocal theories with theories mod- 
eling memory effects according to their basic structural ingredient of nonlocality 
(in space or time): 
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Space 


Time 


Strongly nonlocal 


space integrals 


memory functionals 


Weakly nonlocal 


gradient dependent con- 
stitutive functions 


rate dependent constitu- 
tive relations 


Relocalized 


??? 


internal variables 



The table is self explanatory, what we should observe is, that all of the weakly 
nonlocal approaches can be found in the third row of the second column, that is 
they introduce gradient dependent constitutive quantities (thermodynamic poten- 
tial, entropy current, conductivity tensor, etc..) to generate the nonlocal effects. 
It is conspicuous the lacking of the nonlocal counterpart of internal variables. 

The main purpose of this paper is to confront the variational approach with 
the thermodynamic requirements on the example of Ginzburg-Landau equation. 
In the following section we will review the traditional variational derivation of the 
Ginzburg-Landau equation. In the second section a Ginzburg-Landau like equation 
is constructed from thermodynamic principles preserving the compatibility with 
classical continuum theories. The derivation generalizes the internal variables to 
nonlocal phenomena and based on a general entropy current. With this method the 
and the nonlocal constitutive functions can be generated with the help of a force- 
current system like in classical irreversible thermodynamics . In this way we can 
eliminate several ad-hoc assumptions of the previously mentioned thermodynamic 
approaches. The nonlocality of internal variables is generated by thermodynamic 
requirements and the structure of the theory. In the fourth section a more rigorous 
derivation, based on the Liu procedure, shows the applicability conditions of the 
irreversible thermodynamic treatment. From that considerations we can recognize 
that the previous heuristic thermodynamic derivation is almost completely general 
under some straightforward physical requirements. 



2. Variational derivation of the Ginzburg-Landau equation 

From a thermodynamic point of view the Ginzburg-Landau equation seems to be 
the first nonlocal extension of an evolution equation of an internal variable (order 
parameter, dynamic degree of freedom, etc.). Its traditional derivation is based 
on the introduction of a gradient dependent thermodynamic potential, a so called 
(Helmholtz) free energy functional. In the most simple case this thermodynamic 
potential can be written in the following form 

(1) J {.fio + ^mf)dv. 

where ^ denotes the internal variable, and / is the 'equilibrium part' of the free 
energy. According to intuitive thermodynamic requirements we may assume that 
the equilibrium of the related physical system is characterized by the extremum 
(minimum) of the above free energy functional. In this case we perform a vari- 
ation and get the 'functional derivative' of as a generalized intensive quantity 
conjugated to the internal variable 
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Here the dash denotes a derivative. We can get the stationary Ginzburg- Landau 
equation assuming that this funetional derivative is zero 

no - AAe - 0. 

Introdueing a 'relaxational dynamics' we can generate the time dependent Ginz- 
burg-Landau equation as 

(2) i - -Ik^ = -lif'iO - AAC). 

where I is a positive scalar coefficient (e.g. in case of scalar internal variable and 
isotropic material). We may introduce more general free energies, internal variables 
with different tensorial order, nonisotropic material, etc but the essence of the 
variational derivations remains the same. 

As it was pointed out by Gurtin, the most important problem with that deriva- 
tions is, that they have nothing to do with the balances of the fundamental physical 
quantities. He emphasizes the importance of the separation of balances from the 
constitutive properties: "My view is that while derivations of the form ... are use- 
ful and important, they should not be regarded as basic, rather as precursors of 
more complete theories. While variational derivations often point the way toward 
a correct statement of basic laws, to me such derivations obscure the fundamental 
nature of balance laws in any general framework that includes dissipation." 

There are three specific problems with the variational approach as a result of 
the incompatibility with the structure of continuum physics. 

- The above derivation is nothing to do with a nonnegative entropy production, 
e.g. the sign of I is fixed according to some stability requirements, not by pure 
thermodynamic conditions. 

- The whole equation is not derived from the variational principle, the relax- 
ational dynamics is an additional, independent requirement. 

- The whole procedure restrict how the rate terms can appear in the equation. 
Essentially they can be added by hand waving physical arguments in each 
specific theory, but experiments and theoretical considerations show well, that 
they have a typical form, that is independent of the specific theory and a 
generalized Ginzburg-Landau equation can be written as (see and the 
references therein) 

(3) i - -/(/'(O - AAO + kAi 
Here fc is a positive coefficient. 

3. Thermodynamic derivation of the Ginzburg-Landau equation 

In this case our task is to find an evolution equation of an internal variable that 
corresponds to the requirement of nonnegative entropy production. We suppose 
here that the entropy function depends only on the internal variable and we will 
use the notation := £'s(C)- If C = then r^{0) = 0, because ^ is an internal, 
dynamic variable. As regards the entropy current we apply the previous physical 
assumptions also in this case: if ^ was zero then there is no entropy flow. Moreover, 
in the light of this assumption and according to the mean value theorem, the entropy 
current can be written as a linear function of the derivative of the entropy, as in 
classical irreversible thermodynamics. The coefficient can depend on the internal 
variable, too: 

j,(e)-A(or4. 
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That form of the entropy eurrent was suggested by Nyiri |15 . Therefore the 
entropy production foUows as 



Gs = m + V ■ j. = T^{i + V ■ A) + A ■ vr^ > 0. 

We can recognize a force-current structure. In isotropic materials the two terms 
do not couple and the corresponding Onsagerian equations in the linear approxi- 
mation are 

(4) e + V-A = 

(5) A = 

Eliminating A from (jj) and we get 

(6) i^hT^~V -{hVT^). 

Here we have got an equation that is similar to the Ginzburg-Landau equation 
. However, there are some differences. 

- At the second term of the right hand side, under the space derivatives there 
is instead of ^. However, is a homogeneous linear function of ^, being 
an internal variable. 

- The sign of the material coefHcients is determined by the Second Law, not by 
loose stability considerations. 

- There is no additional rate term at this level of approximation. 

- To extend the derivation to nonlinear and anisotropic cases is straightforward. 
However, the nonlinearities and anisotropics show a different structure, than 
in the original equation. 

This thermodynamic Ginzburg-Landau equation was derived also by Verhas un- 
der some slightly different assumptions, as a governing equation for the transport 
of dynamic degrees of freedom |l^ . 

3.1. Generalized thermodynamic Ginzburg-Landau equation. We can get 

the generalized form of the Ginzburg-Landau equation at the next level of approxi- 
mation. A similar equation was received by Gurtin with the principle of microforce 
balance. In the thermodynamic derivation we consider the previously introduced 
current intensity factor A as an internal variable and follow a similar procedure. 
For the sake of simplicity we suppose, that the entropy function does not depend 
on A, that is we arc not interested in the associated memory effects, we are inves- 
tigating only the nonlocal extension. The form of the entropy current is similar to 
that of the Cahn-Hilliard equation p^ : 

j,(A,e)-AF5 + B(A,e).A. 

Here B is a second order tensor. This form is completely general under the 
conditions of the mean value theorem, if we exploit that there is no entropy flow 
when the internal variable A is zero. Now the entropy production will be 

o-s = r^(C + V ■ A) + A ■ (VF^ + V • B) + B : VA > 0. 

It is straightforward to put down the Onsagerian conductivity equations, but 
after the elimination of B, one cannot simplify them further. Therefore, we will 
treat here only the simplest situation, when the material is isotropic and the ap- 
proximation is strictly linear (the conductivity coefflcicnts arc constants). Now the 
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conductivity equations are reduced to the following form 

(7) e + V-A = hT^, 

(8) A = /2(vre + v-B), 

(9) B = /iVoA + /2(Vo A)* +Z|V- A, 

where /i, ^2, ^3, ^3, are positive, scalar, constant coefBcicnts. Now a simple calcu- 
lation eliminates A and B from the above equations and we get 

(10) i^hT^~l2A{l + hh)r^ + l3Ai 

where I3 ^ + + l^. The last term, that is additional to (H) corresponds to the 
generalized Ginzburg-Landau equation (|^). The positivity (positive definiteness 
in a more general situation) of the material cocfhcicnts is ensured by the Second 
Law. However, we should observe, that this generalization has not changed the 
the characteristic thermodynamic term which differs form the original Ginzburg- 
Landau form: stands under the Laplacian instead of ^. 

We can continue the introduction of new nonlocal internal variables, putting 
B into the basic state space. In this case B becomes internal variable and we 
can introduce a corresponding current intensity factor. Continuing this procedure, 
we can develop a whole phenomenological hierarchy of weakly nonlocal transport 
equations of higher and higher orders. The further research in this direction has a 
special importance for the kinetic theories, because we do not have a well established 
approximation scheme for nonlocal phenomena like the momentum series expansion 
in case of memory effects. The outlined phenomenological hierarchy of nonlocal 
equations can suggest a similar approach for the kinetic equations including the 
sensitive question regarding the closure of the corresponding relations jl^, |l9j. It 
is straightforward to extend the above treatment considering memory and nonlocal 
effects together. 

4. A MORE EXACT DERIVATION OF THE GiNZBURG-LANDAU EQUATION 

There are some basic problems in the heuristic approach of irreversible thermo- 
dynamics that should be addressed in a modern treatment. We can avoid them 
applying a more rigorous form and exploitation of the Second Law (nonnegative 
entropy production) that was applied in the previous section. In this section we 
make a distinction between the state variables and the constitutive quantities at 
the beginning and we apply the Liu procedure 

As we have seen in the heuristic treatment, the thermodynamic Ginzburg-Landau 
equation was received as a general evolution equation for an arbitrary internal 
variable with the requirement of the compatibility with the Second Law. Therefore, 
let us denote our basic state space spanned by an internal variable ^ with Z^. We 
are to find an evolution equation of the internal variable in the form: 

(11) dti + T = 0, 

with the requirement of a nonnegative entropy production 

dtS + V -is > 0. 

Here is the entropy current and dt denotes the partial time derivative. With 
the different notation of the time derivatives we emphasize that in case of moving 
media some further considerations are necessary. The constitution space, where 
the constitutive quantities are defined is spanned by ^ and its first and second 
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gradients (C, V^, V^^). Therefore d = x Lin{Z^,'R.^) x Bilin{Z^,M.^) is the 
constitution space of the nonlocal dynamics of an arbitrary scalar internal variable. 
The constitutive quantities are the entropy, the entropy current and the form of 
the evolution equation (sj^,^). Therefore in the Liu procedure the nonnegative 
entropy production is supplemented by 

dti + J" = 0, 

According to Liu's theorem there exist a F, to be determined from the Liu 
equations, which can be written in a particularly simple form 

dis-T = 0, 

d2S = 0, 

dzs = 0, 

fsjs = 0. 

The dissipation inequality in our case is 

(12) dus ■ Ve + d23s ■■ V^e ^TT>0. 

The solution of the Liu equations gives that s depends only on the internal 
variable ^, F = = -Ds(^) is the derivative of the entropy and js does not depend 
on the second gradient of ^. Unfortunately these considerations do not simplify the 
entropy inequality at all, we may look additional conditions. We write the entropy 
current in the following form 

(13) j,(e,vo = A(e,ve)F(e). 

Considering the Liu equations we can sec again, that this form of the entropy 
current is completely general. Under some differentiability conditions the entropy 
inequality turns out to have the same form that we received with the heuristic 
considerations of the previous section 

A • VF+ (V • A + J^^)F > 0. 

Let us remark, that in (|l^) there arc two constitutive quantities (the entropy current 
and T) and three additive terms. To simplify the inequality for constructing a force 
current system there is no other choice that we have done here: we should unite 
two of the terms with some reasonable physical assumption. Assumption (^|) is 
almost purely mathematical, the physical condition is hidden in the differentiability 
of A. But that is necessary for the construction of constitutive functions, because, 
if A is continuously differentiable we can solve the resulted inequality and get the 
Onsagerian structure. 

Remark 4.1. In this procedure we applied the Coleman- Mizel form of the Second 
Law, with the requirement that the nonnegative entropy production is a consequence 
of pure material properties, and valid for all (continuously) solutions of the evolu- 



tion equation on the internal variable |22 



Let us observe that one of the consequences is that a gradient dependent ther- 
modynamic potential (entropy of free energy) is not compatible with a nonnegative 
entropy production and a relaxational (or any) kind of evolution equation for the 
internal variable. 
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5. Discussion and further remarks 



The applied thermodynamic procedure can be apphed to receive several other 
classical weakly nonlocal equations of continuum physics, e.g. the Guyer-Krumhansl 
equation for nonlocal heat conduction of the Cahn-Hilliard equation for first order 
phase transitions p3| . The Liu procedure shows well the conditions (that are mostly 
mathematical) . 

The differences between the traditional and the thermodynamic form makes pos- 
sible to compare the consequences of the equations experimentally. Understanding 
the (space-time) structure generating properties of the original Ginzburg-Landau 
equation from a thermodynamic point of view seem to be really interesting. Fur- 
thermore, we have seen that the experimentally observed additional rate dependent 
term of the generalized Ginzburg-Landau equation is a natural consequence of the 
thermodynamic approach. More thorough considerations on the possibilities of 
variational principles show well that it would be hard to force that rate dependent 
terms into a variational principle (see e.g. ^). 
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